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Abstract. This talk will discuss how to compute initial quantites in heavy ion colli- 
sions at RHIC (200 AGeV) and at LHC (5500 AGeV) using perturbative QCD (pQCD) 
by including the next-to-leading order (NLO) corrections and a dynamical determina- 
tion of the dominant physical scale. The initial numbers are converted into final ones 
by assuming kinetic thermalization and adiabatic expansion. 



The whole heavy ion physics community has entered into exciting and intense 
period as the first resuhs from RHIC are expected to appear, and the era of LHC 
seems no more very distant. However, despite the sohd status of QCD as the 
theory of strongly interacting matter, many uncertainties remain in the predictions 
for the outcome of the current experiments and various methods have been applied. 
Further experimental data will (hopefully) single out the best candidates for the 
correct approach. 

As the collision energy is increased from that of the SPS, larger intrinsic scales 
are generated and the applicability of perturbative QCD (pQCD) becomes possible. 
The initial particle production is expected to be dominated by minijets, i.e. partons 
with pt ~ 1 . . . 2GeV:^ ^qcd [!]• By assuming independent multiple semi-hard 
parton-parton collisions the average energy carried by minijets with pj- > po at the 
rapidity interval AY in a central b = AA collision is given in leading-order (LO) 
by [2] 
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where T^A(b) is the standard nuclear overlap function and po is the smallest 
transverse momentum scale to be considered. CoUinear factorization is assumed to 
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hold and any effects beyond it are neglected. This perturbative minijet approach 
suffers from two major sources of uncertainties: 

(i) The next-to-leading order (NLO) corrections, which have not been known prior 
to [3,4], but rather have been simulated by ad hoc X-factors. 

(ii) The determination of which value of po to use at, say, RHIC or LHC energies. 
Should one stick to some constant universal value of po or will this parameter 
possess some nontrivial ^/s- and A-dependence? 

In the following we will provide answers to both of these questions and combine 
them to obtain numerical estimates of average transverse energies and charged 
particle multiplicities at RHIC and LHC energies. 

As is evident from formula (1), to deal with the first uncertainty, the relevant 
quantity to compute in NLO pQCD is the cr(£'r), the first moment of the pertur- 
bative Et distribution in a pp-coUision. The infrared safe NLO computation of this 
quantity has been presented in [4] , where the computation was formulated via the 
subtraction algorithm of S. Ellis, Z. Kunszt and D. Soper [5]. 

The Et in central rapidity region is defined to be the total pt entering this 
region and originating from hard subprocessess in which at least an amount of 2po 
of transverse momentum is released. The numerical results are shown in fig. 1(a). 
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FIGURE 1. (a) The NLO a{ET). The LO' stands for the leading order result evaluated 
with 2-loop as and NLO parton distributions, whereas the LO stands for the leading 
order result evaluated with 1-loop and LO parton distributions. The rapidity interval 
is chosen to be the central unit, and the parton distributions used are those of GRV-94 
-set [6]. (b) The average number of QCD-quanta produced with pT > Po and \y\ < 0.5 
as a function of pQ. The saturation scale psat is determined by the points of intersection 
with the dashed curve {pqR\) labelled "saturation". 

To analyze the implications of these numbers, let us define two X-factors: K — 
(fuH NLO) /LO and K' = (full NLO) /LO', the first of these measuring the deviation 
of the NLO results from the consistent LO calculation and the latter measuring 
rather the relative difference between two subsequent terms of the perturbation 
series. The magnitude of these X-f actors is due to a new kinematical region which 
manifests itself only at NLO. If one rejects this new domain, the resulting X-factors 



would be a factor of two smaller than those obtained from fig. 1(a), but then a 
significant amount of perturbatively calculable Et would be neglected. 

For a detailed description of the calculation and issuses such as the scale choices, 
see ref. [4] and references therein. For the purposes to be considered here, the 
sufficient observation is that we now have control over the magnitude of NLO 
corrections which are stable relative to LO results even at few GeV scales, thus 
signalling the applicability of pQCD in this domain. 

Turning to the uncertainty (ii), then, it is clear that as po is decreased the 
cross section, as well as the uncertainty, grows. At certain value of saturation, 
Po = Psat the system becomes very dense and new physics enters [7]. For large 
nuclei and large collision energies this may happen already in the perturbative 
domain, which we concluded on the basis of the NLO analysis to include also the 
few GeV region. Then the corresponding vahies for the number of particles as well 
as for the amount of Et are easily produced via the perturbative computation at 
Po — Psat, which effectively accounts for the contributions of all scales, since the 
partons with px ^ Psat ^-re rare and those with px <^ Psat, altough numerous, 
contribute negligibly to total Et- 

Various ways to determine the actual magnitude of Psat can be conjectured. A 
simple geometric criterion has been presented in [8]. This is based on the idea that 
if one assigns an effective area tt/pq to each gluon produced, then at certain value 
of Po the total area of Naa{po: ^/sAY) gluons produced will exceed the effective 
transverse area 7tR\ of the nucleus. Therefore one can iterate the equation Naa ~ 
Psat-^A to determine psat for given A and ^/s, see fig. 1(b). On the basis of the 
NLO analysis, we take here = 2 to account for the NLO corrections and also 
implement nuclear shadowing via the EKS98 parametrization [9] . 

All the initial quantities then computed can well be fitted by a scaling law of a 
type CA^{^/sy. In particular one finds that 



where the particle and energy densities are evaluated at Psat, and the whole pro- 
duction process is then considered to take place at Tq = 1/Psat- 

These, however, are just the initial numbers at 0.2 (0.1) fm/c at RHIC (LHC), 
and the major problem is how to get from these to the ones at later instants and 
to finally arrive at experimentally visible quantities. Let us therefore assume, not 
completely without reason (sec [8]), that the system is initially thermalized in the 
sense that it possesses a correct ratio of energy per particle as far as the dominant 
gluonic particle content is considered, and expands conserving the total entropy 
S 3.6A^j. As the final particles consist dominantly of pions, for which S fa 4A^/, 
we find that Nf — 0.9Ni. As the initial volume is Vi — 7ri?^Ay/psat, formulae (2) 
give 



Psat/GeV = 0.208A°-i28(^) 

ei/(GeV/fm=^) = 0.103A°-5°^(v^) 
n,-fm^ = 0.370A0-383(^) 
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After the conference the very first measurements by PHOBOS collaboration have 
been announced [10]. According to them, the charged particle multiplicity at midra- 
pidity is dN/dr] = 408 ± 12 ± 30 at 56 AGeV and 555 ± 12 ± 35 at 130 AGeV. The 
framework described here gives Nch = 2/3Nf = 370 and 530 per unit r) respectively, 
when taking into account that number of participants was reported to be 330 for 
^/s = 56 AGeV and 343 for ^ = 130 AGeV, and that dN/dy = l.lddN/dr]. 

The final Et in this scenario is obtained by means of hydrodynamics [8,11] as 

Er^NfX [0.39 + 0.061 ln( A^^ /A)] . (4) 

Numerical values then per unit r] are Et — 260 GeV for -y/s = 56 AGeV and 
Et = 390 GeV for a/s = 130 AGeV using the multiplicites computed above and the 
quoted participant numbers. Measurements of these final transverse energies are 
awaited to appear soon. These measurements will then allow us to draw conclusions 
on the issues such as the true degree of thermalization in the system. 
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